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Abstract 

The scavenging ratios for SO^, SO^ ~, and NO-" vary markedly (over 
two orders of magnitude) between precipitation events and are shown to 
possess log-normal distributions. Using a constant value of scavenging ratio 
as an input parameter into a model, as is generally practiced, will lead to 
significant errors in model calculations. 

Scavenging ratios for 50^ show little spatial variation. The spatial 
variability for NO-," is larger. Also, the scavenging ratios for NO-,* 
suggest that in cloud formation of NO^ may be important. 
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L Introduction 

Current long range transport models employing a linear wet and dry 
scavenging scheme assume that wet deposition of pollutants is given by the 
product of a scavenging coefficient A. (s~ ) and the integral over the mixed 
layer depth (h) of the pollutant concentration in the air. If one defines 
scavenging ratio W as the ratio of pollutant concentration in precipitation 
to Its air concentration, then W equals ^^- where b is the precipitation 
rate. This assumes that the pollutant is uniformly distributed in the mixed 
layer. 

Authors differ in their choices of scavenging parameters. Elliassen 
(197S) assumes A. to be constant. Shannon (1980) defines A as a power law 
function of 3 , Ellenton et al, (1983) assume W to be constant. Unless A(or 
W) is derived from fundamental principles of in-cloud chemical 
transformation and mass transfer of pollutants to water droplets in the 
clouds, and the functional dependence on dominant variables is expressed 
expiicity, one has no guidance on evaluating them by the precipitation 
events. On the other hand, physical and chemical processes in the 
atmosphere dictate that A or W should vary markedly between 
precipitation events. 

Since a linear treatment of the scavenging processes leads to an 
exponential decay type solutions to the pollutant species, the variability of 
A (or W) cannot be simply approximated by using an average value for long 
term simulations as is generally practised by the modellers. It is shown in 
Section 3 of this paper that ignoring the variability of A (or \V) may lead to 
serious errors in the final model computations. 
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The purpose of this paper is to present the distribution of W for SO 2, 
SOf, and NOn' as derived fronn data from a monitoring network (APIOS - 
Acidic Precipitaiton in Ontario Study) in Ontario, Canada and to show the 
implications of this variability of W on model calculations. It is shown that 
W possesses a log-normal distribution. Using an average value of W in 
Lagrangian long-range transport models leads to an under-predlction of 
concentration and deposition of sulphur in remote areas- 

2. Network Descriptiian 

As part of the APIOS, Ontario operates two networks monitoring both 
wet and dry deposition. One network samples on a daily basis (Chan et al, 
1982) and the other one 28-day basis (Chan et al. 198^). This paper presents 
data obtained in the former network. 

The daily network has been in operation since July, 1980 at one site 
(Dorset) in central Ontario. It was expanded to 16 sites in 1981 which are 
cluster d into 4 groups with spatial separation of more than 1,000 km. Wet 
deposition is collected at all 16 sites with Aerochem Metrics wet-only 
samplers at all times except during the winters of 1980/81 and 1981/82 
when Z^ hour bulk samplers were used. Dry deposition is estimated by 
measuring ambient air concentrations using a lo-vol filtration technique at 
I station of each of the ^ clusters. Samples collected on the teflon 
prefiiters are analyzed for particulate SO^ ', NO-," and NH^"^ whereas 
those on the downstream nylon and impregnated K2CO^-giycerol Whatman 
^i filters are analyzed for gaseous HNO^ and SO2 respectively. Both daily 
wet and dry samples are collected typically around S a.m. Data presented 
here are obtained from Charleston Lake (44°29'54"N, 76^02'30"W), Dorset 
(^5^13'23"N, 78^55'^9"W), Longwoods (42^53'02"N, 81°28'50"W) and 
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Fernberg (47^36'51"N and 91^29'26"W). The 3 former sites are located in 
Ontario. The last site is located in Minnesota where the Ontario 
instruments are co-located with those operated by the United States 
Environmental Protection Agency at Duluth. 

3. Reajjts and Discussion 

Only stations with both precipitation (mg 1" ) and air (ug m" ) 
concentration results are examined for scavenging characteristics of the S 
and N compounds in 6 different combinations: 



2_ 
W I r S-SO^^ " (precip. 



'^ 



S-SO^ (air) 



2 

W2 = S - SO^ ~ (precip.) 



S - SO2 (air) 



W^ = S - %0^' (precip.) 



2- 
S - SO^ " (air) + S - SO2 (^i"^^ 



W^ = N - NO3" (precip.) 
N - MO-j" (air) 



W^ = N - NO3" (precip.) 
N - HNO3 ^^^^^ 



W^ - N - NO3' (precip,) 



N - NO3' (air) + N - HNO3 (air) 



- ^ - 

2- - . 

The above combinations are necessary as SO^ and NO^ in 

2- - . 

precipitation do not originate entirely from particulate SO^ and NO^ in 

air. Some unknown portions originate also from SO- and HNO-, (and 
possibly NOj^). W^ W2, W^ and W^ yield upper bounds of scavenging ratios 
for SO.^", SO2, NO^" and HNO^. In order to account for the fraction of 
precipitation concentration due to other co-existing species in the air, W 
shouid be scaled accordingly. In the case of SO^" and SO2 in air, Henry's 
law calculations could be appOed (Barriej 1981). It is recommended that 
factors 0.9 and 0.1 should be applied to W. and W2 as a first order 
approximation as it has been found that SO2 correction is less than 10% 
(Barrie, 1983), The correction for the N species is not well understood, 
however, Huebert et al. (1983) have reported a ratio of W^ to W, being five 
by snow scavenging. 

The calculated ratios are tested for normal and log-normal 
distributions using Kolmogorov-Smirnov statistics (Pollard, 1977), when the 
number of data points is greater than 50 (otherwise the Shapiro-Wilk 
statistics is used). Basically, this test compares the largest absolute 
cumulative difference between the tested sample distribution and the 
normal distribution with the Kolmogorov-Smirnov^ statistics (or Shapiro- 
Wilk (w) statistics) of the same sample size. When the test value (D ^ J is 

Stat 

greater than the critical value (DQjg) at 5% significance level, it indicates 
that the sample possibly does not come from a normal distribution. In other 
words, we may conclude that the sample distribution is not normal with a 
5% chance of making a wrong conclusion. However, only when the test 
value (D ^) is well below the critical value (Dq q J, we may have high 
confidence in concluding that the sample distribution is normal. Results of 
the test (Tables 1 to 6) indicate that the observed scavenging ratios follow 
a log normal distribution well. The following points are noteworthy. 
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a» There is little spatial variability In SO^ scavenging ratio. 



h. The spatial variation in NO-" scavenging ratio is larger reflecting 
non-uniform change in concentration of NO^" in precipitation and 
NO-' and HNO3' in air geographically. 

c. Using Junge's scheme (1963) of rainout efficiency €., or W. x 10" iin 
mir notation) Hidy (19S2) indicated that rainout efficiency could be 
identified with the scavenging ratio based on ground level 
precipitation and air measurements. With simplifying assumptions 
including (i) below cloud scavenging is small compared with in-cloud 
scavenging processes; (ii) clouds have an average 1 g m' liquid-water 
content; (iii) precipitation composition is identified with cloud water 
composition and can be related to ambient air concentration which is 
taken to be constant with height, C- is expected to fall in the range 
of zero to unity. The fact that W, (x 10" ) is greater than unity most 
of the time suggests that NO^~ in precipitation cannot be simply 
explained by its air concentration and that in-cioud oxidation 
formation of NO-," may be important. 

4. Implication and Appiication to McxJels 

Models using a linear scavenging scheme are typically applied for 
deposition calculation over a period of one month to one year. Let us 
suppose that over this time period n. trajectory end points arrive at a 
receptor during precipitation events after a travel period of kAt hr, where 
At is the time step of computation in the model and k is the number of 
time steps. 
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The n, trajectory end points might encounter several precipitation 
events over this time period. The question is, what is the effect on model 
calculation of wet deposition at the receptor by using a constant versus a 
variable W where W is allowed to vary between the precipitation events? 

To answer this question we note that the mean deposition rate of the 
ensemble, obtained by averaging over n, trajectories, is given by the 
followingi 



D=[lr^2^. ^ ^4^4!^ 



(1) 



2- 
where D is the average rate of wet SO^ deposition, C- the mean SO2 air 

concentration, C^ the mean SO^ ' air concentration, A. 2 the mean 

scavenging coefficient for SO^f ^ [, the mean scavenging coefficient for 

2- 
SO^ and h the mixed layer height. 



We are concerned here only about the deposition at a fixed receptor 
averaged over several events. This is different from the average effects 
along a single trajectory. 

The critical parameters In (1) are C^ and C^ which depend on how 
A and A vary between the trajectories encountering precipitation 
events. 

C2 and C^ are given as follows: 






C2= I Ci(./v^T;pcA,}c(A^ (2) 
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*= j j^/^2^:; A^+)|?(A^,A^)4A^dA^ 



and C = 

o o 



(3) 



Here Co and C^ are expressed as functions of A 2 and A^. We assume that 
A and A. are determined by the same physical processes and therefore, 
A-/A. = o( , a constant. The probability density functions pCA^) and p(A^) 
possess similar shapes differing only by the constant factor oC . 

It is general practice in modelling to approximate C2 and C^ by the 
following: 



C2 = C2 CA^t) 



(^) 



and C ,. = C^ {A^Js.X) ^ ^^^ 



'% 



To illustrate the difference between (2) and (^) as well as (3) and (5), 
let us consider <Z^. Its ratio R as calculated by (2) and (^) is given by the 
following: 

(qCA,+ )^CA,)dA^ 

R= —1 (6) 

It is noted that since A-2 ranges from to o^ and C2 (-A-t) is a 
decaying function of A2t, there is a range of values of A 2 for which C^ 
{K^ decays slower than Z^ (A^t). Thus the ratio, CJ.K^)lC2 (^'^)j will be 
an incresing function of t in the range, ^ A S Jl . In other words, R 
win increase with t. 

This implies that we will underpredict C- by using equation (^), as 
travel time increases. Similar conclusions can be reached for C^- In areas 



considerably removed from source areas we will underpredict the 
concentration of SO^ by using an average value of -A.^ in modelling rather 
than a variable A-. Precisely where we will begin to observe this 
difference and its magnitude will depend on the value of A- and the 
frequency with which precipitation occurs along the trajectories. This will 
be further investigated in a Lagrangian long range transport model. 

^ Conciusions 

The following conclusions are drawn from this study. 

L The scavenging ratios of S and N species are highly variable 

with a range spanning two orders of magnitude. 

:1» The scavenging ratios are log-normally distributed. The 

reason for this distribution is as yet unclear. 

3» Using an average value of the scavenging ratio in linear long 

range transport models may lead to serious errors in model 
predicuons at remote areas. The use of a variable scavenging 
ratio based on the distribution presented here is 

recommendedt 

2- 
i* SO^ scavenging ratios show little spatial variability. NO^ 

scavenging ratios show greater spatial variability. 

S, The larger the^n unity (W^ x 10' ) scavenging ratios for NO," 

suggest that in-cloud formation of NO^' nnay be important. 
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Table 1 

2- 2- 

Sunmary of Scavenging Ratio (W, = S-SO^ (precip. )/S-SO^ (air)) Results 

Non Transformed Log Transformed W-(x 10 ) 

Data Data 

Station N D.f,* 0« nc 0^.^. Dn nn Geom. Mean Geom. Std. Dev 

Stat 0.05 Stat 0.05 

Charleston Lake 216 0.335 0.0924 0.0526 0.0924 0.86 3.02 

Dorset 330 0.385 0.0748 0.0732 0.0748 0.88 3.38 

Fernberg, Mn 93 0.379 0.1408 0.0757 0.1408 0.70 3.20 

(N - 92) 

Longwoods 197 0.481 0.0968 0.0700 0.0968 0.82 2.86 



Table 2 
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Summary of Scavenging Ratio (W^, = S-SO^ (precip. )/S-SO^ (air)) Results 



Non Transformed 
Data 



Log Transformed 
Data 



W^Cx TO"®) 



•Station 



Charleston Lake 



n. 

213 



Stat 



0.385 



%.05 ^stat 



0-0931 0.0417 



'0.05 



0.0931 



Geom. Mean Geow, Std. Dev 



0.46 



4.63 



Dorset 



Fernberg, Mn 



Longwoods 



284 



71 



189 



0.413 0.0806 0.0361 



0.349 0.161 0.0704 



0.336 0,0988 0.0443 



0.0806 



0.161 



0.0988 



0.50 



0.76 



0.34 



4.80 
4.70 
3.91 



Table 3 

Summary of Scavenging Ratio (W^ = S-SO^ (precip. )/(S-S02 (air) + S-SO^ (air))) Results 

Non Transformed Log Transformed W-(x 10~ ) 

Data Data 

Station N 0^.,,. D^ nc D^i-,4- On nn Geom. Mean Geom. Std. Dev 

Stat 0.U5 Stat U.U5 

Charleston Lake 216 0.254 0.0924 0.067 0.0924 0.26 3.36 

Dorset 294 0.372 0.0792 0.0465 0.0792 0.29 3.71 

Fernberg, Mn 93 0.415 0.1408 0.086 0.1408 0.37 3.56 

(N = 92) 

Longwoods 190 0.275 0.0985 0.0587 0.0985 0.22 3.19 



Table 4 



Summary of Scavenging Ratio (W^ = N-NO^ (precip.pt)/N-NOo (air)) Results 



Station 
Charleston Lake 



n 

200 



Non Transformed 
Data 



stat 



0.363 



0.05 



0.0960 



Log Transformed 
Data 



""stat 

0.0442 
(N = 199) 



'0.05 



0.0962 



W,(x 10-®) 



Geom. Mean Geom. Std. Dev 



5.64 



4.23 



Dorset 



263 



0.365 0.0837 



0.0450 
(N = 261) 



0.0841 



18,63 



5.71 



Fernberg, Mn 



Longwoods 



45 



192 



W* = 0.600 



* The Shapiro-Wilk Statistics 



W - 0.965 
(N = 44) 



0.425 0*0980 0.0843 



0.0980 



7.61 



2,41 



5.46 



4.28 



Table 5 



Summary of Scavenging Ratio (Wj- = N-NOo" (precip. J/N-HNO^, (air)) Results 



Non Transformed 

Data 



Log Transformed 
Data 



W^(x 10"^) 



Station 



Charleston Lake 



N 
212 



stat 



0.246 



^0.05 ^stat 



0.0933 0.0614 



0.05 



0.0933 



Geom. Mean Geom. Std. Dev 



1.49 



2.58 



Dorset 



335 



0.429 0.0742 



0.0751 0.0745 
(N = 332) 



1.82 



2.98 



Fernberg, Mn 



92 



0.382 0.142 



0.0890 0.144 
(N = 89) 



2.51 



3.40 



Longwoods 



193 



0.302 



0.098 0.0603 



0.098 



1.38 



3.09 



Table 6 



Summary of Scavenging Ratio (Wg = N-NO^ (precip. )/(N-NO^ (air) + N-HNO^" (air))) Results 

Non Transformed Log Transformed Wc(x 10" ) 



Data Data 



6 



Station N D^^^^ D^ ^^ D^^^^ Dp ^^ Geom. Mean Geom. Std. Dev 

Charleston Lake 213 0.256 0.0931 0.0600 0.0931 1.06 2.55 

Dorset 332 0.395 0.0745 0.0652 0.0745 1.53 2.86 

Fernberg, Mn 88 0.343 0.1448 0.0822 0.1448 2.09 3.42 

Longwoods 198 0.301 .0965 0.0669 .0965 0.74 2.67 
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